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Abstract: We present the double real radiation contributions to the ti hadronic produc- 
tion cross section stemming from the partonic process gg — t- tiqq. We explicitly construct 
the antenna subtraction terms for this gluon-gluon initiated process starting from the dou- 
ble soft behaviour of the double real radiation matrix elements using soft currents. Those 
subtraction terms, given in leading and subleading colour contributions, require the use of 
new genuine NNLO four-parton antenna functions involving massive fermions. Those are 
also presented together with their infrared limits in this paper. We checked the validity 
of our subtraction terms numerically by showing that the ratio between the real radia- 
tion matrix elements and the subtraction terms approaches unity in all single and double 
unresolved regions of phase space. 

Keywords: QCD, Jets, Colhder Physics, NLO and NNLO calculations with massive 
particles. 



Contents 



1. Introduction 2 

2. Antenna subtraction for double real radiation with massive final state 
fermions 5 

2.1 Double real radiation eontributions to heavy quark pair production 6 

2.2 Single unresolved radiation subtraction: dtr^^^^ 8 

2.3 Colour-connected double unresolved radiation subtraction: dcr^^^Q 9 

2.4 Colour-unconnected unresolved radiation subtraction: da^^^^Q 10 

3. Antenna functions 11 

4. Infrared limits 16 

4.1 Single unresolved factors 16 

4.2 Double unresolved factors 17 

4.2.1 Colour-unconnected pairs of unresolved particles 18 

4.2.2 Colour-connected pairs of unresolved particles 18 

4.3 Single and double unresolved limits of four-parton antennae 21 

5. Top quark pair production at the LHC 22 

5.1 Conventions 22 

5.2 Double real radiation contributions 23 

5.3 Subtraction terms 26 
5.3.1 Double soft behaviour of amplitudes 26 

S h 

5.4 The construction of da ' ^/=, _ 28 

gg^QQqq 

5.5 Construction of the do"'^''^ „ =, _ and dcr'^''^ _ subtraction terms 30 

gg^QQqq gg^QQqq 

6. Numerical results 32 

6.1 Double soft limit 33 

6.2 Triple collinear limit 33 

6.3 Double collinear limits 33 

6.4 Final-final single collinear limit 34 

6.5 Initial-final single collinear limits 35 

7. Conclusions 35 

8. Acknowledgements 36 



- 1 - 



A. 



Three-parton antennae 

A.l Massive final-final antennae 
A. 2 Massless initial-final antennae 
A. 3 Massive initial-final antennae 
A.4 Massless initial-initial antennae 



36 



37 
37 
38 
40 



B. 



Four-parton antennae 

B.l Massive final-final B-type antennae 

B.2 Massless initial-initial G-type antenna 



40 



41 



41 



1. Introduction 

With a mass rrit = 173 it 1.3 GeV, the top quark is the heaviest quark produced at colliders 
and due to its very large mass it decays before it hadronises. Top quarks [1,2] are measured 
through their decay into a bottom quark and a subsequently decaying W^-boson, yielding up 
to six-jet final states for top quark pair production. Meaningful searches for these signals 
require not only a very good anticipation of the expected signal, but also of all Standard 
Model backgrounds yielding identical final state signatures. 

By studying the properties of top-quarks in detail, it is hoped to eludicate the origin 
of particle masses and the mechanism of electroweak symmetry breaking. Since its dis- 
covery at the Fcrmilab Tevatron, a number of its properties (mass, couplings) have been 
determined to an accuracy of ten to twenty per cent. With the large number of top quark 
pairs expected to be produced at the LEG [3], the study of its properties will become 
precision physics. In particular the tt production cross section is expected to be measured 
with an accuracy of the order of five percent. These very precise measurements have to be 
matched unto equally accurate theoretical predictions. Fixed-order calculations of these 
observables need to be considered at least at the next-to-leading (NLO), if not even at 
next-to-next-to-leading order (NNLO) in perturbative QCD. 

For any hadronic observable, a theoretical prediction is obtained at a given order in 
as when all partonic channels contributing at that order to the partonic cross section are 
summed and convoluted with the appropriate parton distribution functions. Beyond lead- 
ing order, these partonic channels contain both ultraviolet and infrared (soft and collinear) 
singularities. The ultraviolet poles are removed by renormalisation, while collinear poles 
originating from the radiation of intial-state partons are cancelled by mass factorisation 
counterterms. The remainining soft and collinear infrared poles cancel among each other 
only when all partonic channels are summed over [4, 5] . 

NLO predictions for top quark pair production cross sections have been known already 
for some time [6-8]. Next-to- leading-logarithmic resummation (NLL) are also known [9-11], 
and even the NNLL resummation effects have been completed in [12]. These predictions 
lead to a theoretical uncertainty of the order of ten per cent. The same precision is available 
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for single top quark production [13], top-pair-plus-jets production [14-16] and for top-pair- 
plus-bottom-pair production [17,18]. 

The top quark appears as virtual particle in hadron collider processes and due to the 
small ratio between the top quark width and its mass, it is possible to factor the cross 
section of processes involving top quarks into the product of the production cross section 
for on-shell top quarks and the top quark decay width. Most of the calculations mentioned 
above are performed for on-shell top quark pair production and wc shall also follow this 
approach in this paper. Only most recently, the decay of the top quark has been included 
in NLO calculations [6, 7, 19, 20], leading to a similar theoretical accuracy of the NLO 
predictions. 

At NNLO, intermediate results concerning the top-quark pair production cross sec- 
tion predictions have become available recently: Most notably, the inclusive total hadronic 
tt production cross section induced by the all-fermion partonic processes has been com- 
puted [21,22]. NNLO calculations involving massive quarks require the same ingredients as 
their massless counterparts. Three classes of contributions enter: double real dcr^^, mixed 
real- virtual, dcr^^, and double virtual contributions da^^ . 

Recent progresses has been accomplished concerning the two-loop contributions and 
many parts of the full top-pair NNLO matrix elements are known. Part of these two- 
loop virtual corrections arc built with products of one-loop virtual amplitudes. Those 
corrections have been computed in [23-25]. Concerning the two- loop virtual corrections 
built with product of two-loop and tree-level amplitudes, the situation is different. The 
two-loop virtual corrections for the processes qq tt and gg tt are not fully available at 
present. A purely numerical evaluation of the quark-initiated process [26] could be partly 
confirmed by analytical results [27-29], which were most recently extended also to the 
gluon-induced subprocess. The infrared structure of mixed-real virtual contributions has 
been studied in [30]. 

While infrared singularities from purely virtual corrections are obtained immediately 
after integration over the loop momenta, their extraction is more involved for real emission 
(or mixed real-virtual) contributions. There, the infrared singularities only become explicit 
after integrating the matrix elements over the phase space appropriate to the differential 
cross section under consideration. Since hadronic observables depend in general in a non 
trivial manner on the experimental criteria used to define them, they can only be calculated 
numerically. The computation of hadronic observables including higher order corrections 
therefore requires a systematic procedure to cancel infrared singularities among different 
partonic channels before any numerical computation of the observable can be performed. 

Subtraction methods explicitly constructing infrared subtraction terms which coincide 
with the full matrix element in the unresolved limits and are sufficiently simple to be 
integrated analytically are well-known solutions to this problem. Starting from methods 
for subtraction at NLO [31-35], several NNLO subtraction methods have been proposed 
in the literature [30,36 48], and have been worked out to a varying level of sophistication. 

A subtraction method based on the transverse momentum structure of the final state 
has been proposed in [39]. Using this formalism, NNLO results were obtained for Higgs 
production and vector boson production [49-51]. This formalism has also recently been 
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applied to compute NNLO corrections to associated WH production [52] and photon pair 
production [53]. 

Other approaches to perform NNLO calculations of hadronic observables, concern the 
use of sector decomposition (possibly in combination with subtraction) . The sector decom- 
position approach is a numerical method which avoids the need for analytical integration, 
and which has been developed for virtual [54-57] and real radiation [57-61] corrections to 
NNLO. It has been applied to several observables like Higgs production [62-65] or vector 
boson production [66] . 

It is the combination of subtraction with sector decomposition [30, 40, 41] which was 
also recently applied in the calculation of NNLO corrections to the all fermion-pair initiated 
processes in top quark pair production [21,22] mentioned above. 

Employing a subtraction method, the NNLO partonic cross section for a hadronic 
observable can be written as [67, 68] 



Jd^m+l Jd^m+1 



^NNLO 

Cm+2 



MF,1 
NNLO 

+ / ^^NNLO+ / ^^NNLO- (1-1) 

In this equation, da^jyLO denotes the subtraction term for the (m + 2)-parton final state 
which behaves like the double real radiation contribution da^^^Q in all singular limits. 
Likewise, da^f^j^Q is the one-loop virtual subtraction term coinciding with the one-loop 
(m + l)-final state da^^j^Q in all singular limits. The two-loop correction to the (m -|- 2)- 
parton final state is denoted by da^^j^Q. In addition, as there are partons in the initial 
state, two mass factorisation contributions, da^^j^'^Q and da^^j^'^Q, for the (m -|- 1)- and 
m-particle final states respectively, need to be taken into account. 

In this paper, in order to evaluate the double real contributions to the hadronic top 
quark pair production arising from the partonic process gg — t- QQqq, we shall employ 
the antenna subtraction formalism generalised to hadronic observables involving massive 
fermions developed in [69]. 

Originally, the antenna subtraction method has been developed for the production 
of massless partons in electron-positron annihilation [67,70,71]. The method was then 
generalised in order to deal with coloured initial states and massless final state partons 
in [68,72-77]. The framework for the construction of NNLO antenna subtraction terms for 
hadronic jet observables has been set up in [68, 75] in the context of a proof-of-principle 
implementation of the contribution of the purely gluonic contributions to di-jet production 
at hadron colliders. 

The formalism has also been expanded to include massive fermions. In [78-80] it was 
extended at the NLO level and employed to compute the real NLO corrections to the tt 
and tt -\- jet hadronic production cross section. In [69] we extended the method to the the 
NNLO level, and derived the double real emission corrections to tt hadronic production 
cross section coming from all partonic channels involving fermions only. The aim of this 
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paper is to show how, within the framework presented in [69], partonic processes involving 
gluons can also be evaluated. 

For QCD observables involving massive final state particles, fewer subtractions terms 
are required since the real matrix elements develop singular behaviours in fewer regions 
of phase space than in the massless case. Indeed, for those observables, QCD radiation 
from massive particles can lead to soft divergencies but cannot lead to strict collinear 
divergencies, since the mass acts as an infrared regulator. As a consequence, less divergent 
contributions arise and the infrared structure of such processes is expected to be simpler. 
The kinematics is however more involved due to the finite value of the parton masses and 
the integration of the subtraction terms is more difficult. In addition, in calculations of such 
observables, logarithms involving the ratio of the scales present in a given reaction occur. 
These finite logarithms are related to the singular behaviour of the matrix elements in the 
massless limit [78,80-82]. Depending on the kinematics of the reactions involved, these 
logarithms, although finite, may or may not be enhanced. For the kinematical situation 
under consideration in this paper, i.e. the production of a heavy fermion pair whose mass 
is of the the same order as as the partonic center-of-mass energy, we shall ignore these 
finite logarithms as also done in [21,22], since they are not expected to be large. 

The paper is organised as follows: In section 2 we briefly review the generalisation of 
the antenna subtraction method presented in [69] enabling the construction of subtraction 
terms for the double real radiation contributions to tt production in hadronic collisions. In 
section 3 we give a list of all genuine new NNLO massive four-parton tree-level antennae 
required in our calculation while in section 4 we list the behaviour of these antennae in their 
infrared limits. In section 5, we first present the double soft behaviour of the real matrix 
element associated to the process gg — t- QQqq. This enables us to construct the subtraction 
term capturing the double unresolved features of the double real matrix elements explicitly. 
All subtraction terms required to capture the single and double unresolved behaviour of 
the double real matrix elements in leading and subleading colour will be given in that 
section too. In section 6, we test the validity of the subtraction terms. We check that the 
ratio between the real radiation matrix elements and the corresponding subtraction terms 
approaches unity in all single and double unresolved limits. Finally, section 7 contains our 
conclusions and an outline. Two appendices are also enclosed: Appendix A contains a list 
of all tree-level three-parton antennae together with their infrared limits, while the known 
four-parton antennae that are used in the construction of our subtraction terms are given 
in Appendix B. 

2. Antenna subtraction for double real radiation with massive final state 
fermions 

The extension of the antenna subtraction method for double real emission corrections 
to observables involving heavy fermions in the final state has been extensively discussed 
in [69]. In this paper we consider the double real emission contributions to the hadronic 
tt production stemming from the partonic process gg — )• QQqq. The treatment of this 
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contribution does not require any further extensions of the method. We shall therefore 
give here only a brief summary of the method and refer the reader to [69] for more details. 

Within the antenna subtraction formalism, which uses colour ordering properties of 
amplitudes in a crucial manner, the subtraction terms are constructed using the funda- 
mental factorisation properties of QCD amplitudes and phase spaces in their collinear and 
sof limits. These terms rely on the following essential ingredients: 

(1) A set of antenna functions of various types, which capture all unresolved radiation 
emitted between two hard partons, the radiators. These can be either massless or 
massive and depending on where the two radiators are located, in the initial or in 
the final state, we distinguish three types of antennae: final-final, initial-final and 
initial-initial. 

(2) An exact momentum conserving and Lorentz invariant phase space factorisation 
based on 3 — >^ 2 and 4^2 mappings between on-shell partons in all three con- 
figurations. 

(3) Phase-space mappings defining the momenta present in the factorised form of the 
matrix elements in terms of the original momenta present in the real radiation matrix 
elements for which the subtraction terms are constructed. 

The subtraction terms are then built with products of antenna functions and reduced 
matrix elements squared in all three configurations: final-final, initial-final and initial- 
initial. 

Compared to the massless antenna formalism, the presence of massive partons in the 
final states modifies the subtraction terms in a non-trivial way. Parton masses lead to mod- 
ified kinematics and have to be taken into account for the phase space factorisations [69]. 
Furthermore, non-vanishing masses also modify the soft behaviour of the real radiation 
matrix elements as will be explained in section 5. 

2.1 Double real radiation contributions to heavy quark pair production 

The method developed in [69] enables us to describe the computation of the double real 
radiation corrections to the process 

pp ^ QQ + {m - 2)iets. (2.1) 

To set the notation, the leading order (LO) partonic contribution to this process reads 

daLo{Pl,P2) = MlO XI d^rn{PQ,PQ,Pb, ■ ■ ■ ,Pm-2;Pl,P2) 
m-2 

X-^^\Mm+2{PQ,PQ,P5,- ■ ■ ,Pm-2;Pl,P2)f jjhHPQ^PQ^P5, ■ ■ ■ ,Pm-2), (2.2) 

where pi and p2 are the momenta of the initial state partons, the massive partons Q and 
Q have momenta pQ and pQ, while the momenta of the remaining (m — 2) massless final 
state partons are labelled P5 ■ ■ ■ Pm-2 ■ Sm-2 is a symmetry factor for identical massless 
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partons in the final state. Jrn^\pQ,PQ,P5, ■ ■ ■Pm-2) is tlic jet function. It ensures that out 
of (m — 2) massless partons and a pair of heavy quarks Q and Q present in the final state at 
the partonic level, an observable with a pair of heavy quark jets in association with (m — 2) 
jets is built. At this order, each massless or massive parton forms a jet on its own. The 
normalization factor Mlo includes all QCD-independent factors as well as the dependence 
on the renormalised QCD coupling constant and is process dependent. denotes 
the sum over all configurations with (m — 2) massless partons. is the phase space for 

an m-parton final state containing (m — 2) massless and two massive partons with total 
four-momentum Pi +p^. In d = 4 — 2e space-time dimensions, this phase space takes the 
form: 

X 2^9^ • • • 2^^ S^p.^p.-p,-p,-p,-... . (2.3) 

In eq.(2.2), |A1m+2p denotes a colour-ordered tree-level matrix element squared with 
m final state partons, out of which two are massive and two are initial state partons. 
These terms only account for the leading colour contributions to the squared matrix ele- 
ments, since subleading colour contributions involve in general interferences between sub- 
amplitudes with different colour orderings. However, to keep the notation simpler we 
denote these interference terms also as |7W,n+2p- 

The double real emission contributions to pp — > QQ + {m — 2)jets at the partonic level 

read 

<i^NNLo{Pl^P2) =J^NNLO '^d^m+2{PQ^PQ^P5, ■ ■ ■ , Pm\Pl,P2) 

m 

^-^\Mm+iijPQ,PQ,Pf„ ■ ■ ■ ,Pm;Pl,P2)f J^t^^'^HpQ^PQ^PS, ■ ■ ■ ,Pm)- (2-4) 
= ^NNLO X] d^m+2(P3, • • • ,Pm+i\Pl,P2) 



X^\Mm+i{p^, ■ ■ ■ ,Pm+4, ;^>1,P2)P J^+^^PS, • • • ,Pm+4), (2.5) 

where the last line is obtained by relabelling all final state partons. In this equation, the 
jet function J^'^^^ ensures that out of m massless partons and a QQ pair, an observable 
with a pair of heavy quaxk jets in addition to (m — 2) jets, is built. The next-to-next-to 

leading order normalisation factor A/jy^^Q includes all QCD-independent factors as well 
as the dependence on the renormalised QCD coupling constant a^. It is related to the 
normalisation factor present at leading order, Mlo, which depends on the specific process 
and parton channel under consideration. For the process under consideration in this paper, 
this relation will be specified in section 5. 

This NNLO contribution to the hadronic tt production cross section develops singu- 
larities if one or two final state partons arc unresolved (soft or coUincar). Depending on 
the colour connection between these unresolved partons, the following configurations must 
be distinguished [67,68]: 
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(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour-connected but share a common radiator 
(almost colour- unconnected). 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour-unconnected) . 

(e) Compensation terms for the over-subtraction of large angle soft emission. 

This separation among subtraction contributions according to colour connection is 
valid in all final-final, initial-final or initial-initial configurations and, in any of them, the 
subtraction formulae have a characteristic structure in terms of the required antenna func- 
tions. This antenna structure has been derived for processes involving only massless par- 
tons, for the final-final and initial-final cases in [67, 73] and [68] for the initial-initial case. 
Note that the presence of massive partons in the final state does not modify the general 
structure of the subtraction terms required to match the unresolved features given above 
as explained in [69]. 

For the partonic process gg QQqq considered in this paper, the configuration (c), 
the almost colour-unconnected case, and the configuration (e), regarding the treatment of 
large angle soft radiation, do not occur. On one hand, the colour structure of the amplitude 
for this process does not allow configuration (c) and, on the other hand, the absence of final 

state gluons forbids configuration (e). In the following, we shall therefore restrict ourselves 
to describe how the configurations (a), (6) and (d) are dealt with here while leaving the 
discussion of configurations (c) and (e) to be treated elsewhere. 

2.2 Single unresolved radiation subtraction: da^^^Q 

The general form of the subtraction terms for single unresolved radiation associated to QQ 
production in association with (m — 2)- jets can be taken from the expressions derived 
in [69,78] in all three final-final (f-f), initial-final (i-f) and initial-initial (i-i) configurations. 

These types of subtraction terms, denoted by da^^^^Q, involve the phase space for the 
production of (m -|- 2) partons, with two of them being massive, the three-parton 

antenna functions appropriate to the given configuration X?.^ (f-f), Xfjj^ (i-f) and Xf}; j 
(i-i), the colour-ordered reduced (m -|- 3)-parton amplitude squared |A4m+3p (with one 
parton less than the original amplitude squared) and the jet function J^^^\ Depending 
on the configuration considered, the momenta present in the reduced matrix elements and 
in the jet function arc defined with a particular 3 — )• 2 mapping. These mappings have be 
defined previously in [67-69,72]. 

In all configurations (f-f, i-f and i-i), the products of three-parton antennae and reduced 
matrix elements present in the da^^^Q subtractions terms, coincide with the double real 
radiation matrix element squared given in (2.5), when parton j is unresolved. At NNLO, 
however, the jet function Jm allows one of the (m -|- 1) final state momenta to be 
unresolved and the reduced matrix elements present in the subtraction terms can have a 
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further unresolved parton. In this hmit, the subtraction term d(T^^^Q does not coincide 
with the double real emission contribution given in eq.(2.5) and, as explained in detail 
in [67, 68] , the resulting spurious double unresolved singularities must cancel against pieces 
of the subtraction terms coming from other configurations. In this paper, these other 
configurations will be either (b) or (d), associated to two unresolved particles which are 
colour-connected or colour-unconnected respectively. 

To define the integrated forms of these subtraction terms dcr^^^Q involving only three- 
parton antennae, an appropriate phase space factorisation needs to be considered in all 
three configurations. This factorisation has been derived in the massless case in [72, 74] 
and in the massive case in [78,80]. 

S b 

2.3 Colour-connected double unresolved radiation subtraction: da^Jj^^^Q 

When two unresolved partons j and k are adjacent and colour-connected to two hard 
radiators labelled i and I, we build our subtraction term starting from the four-particle 
tree-level antennae in the corresponding configurations Xij^i (f-f), Xijkl (i-f) and Xujk (i- 
i). By construction, these subtraction terms contain all colour-connected double unresolved 
limits of the (rn + 4)-parton matrix element associated with partons j and k unresolved 
between radiators i and I. However, the four-parton antennae can also be singular in single 
unresolved limits associated with j or k where they do not coincide with limits of the matrix 
elements. To ensure that these subtraction terms are only active in the double unresolved 
limits of the real radiation matrix elements, we subtract the appropriate single unresolved 
limits of the four-particle tree-level antennae using products of two tree-level three-particle 
antennae in the appropriate configurations. 

More precisely, the complete do-^^^Q subtraction terms involve: the phase space for 
the production of (m + 2) partons, with two of them being massive, the colour- 

ordered reduced (m -I- 2)-parton amplitude squared |A^to_|_2|^, (with two partons less than 
the original matrix element squared), the four-parton antennae and products of three- 
parton antennae appropriate to the given configuration. In addition, it also contains they 
jet function jj^^ . The explicit expressions for these subtraction terms in final-final, initial- 
final and initial-initial configurations read, 




^ ^ {^ijkl ~ -^ijk^IKl ~ -^jkl^iJl) 



jk 



,PI,PL,---), 



(2.6) 




1 




i=l,2 jk 



X \M^+2{...J,L,...)\^J^\...,pL,...), 



(2.7) 
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perms m+2 
il=l,2 jk 

X \Mm+2{...J,L,...)\^jt\p3,...,Pm+A), (2.8) 

where the sum runs over all colour-adjacent pairs j, k and implies the appropriate selection 

of hard momenta i,l which, as usual, have three possible assignments of radiators. In all 
cases the (m + 2)-parton matrix clement is evaluated with new on-shell momenta given 
by a momentum mapping appropriate in each configuration. Those 4 — >■ 2 mappings have 
been derived in [68,72,80] in the massless case. In [69] we checked that the same mappings 
can also be applied in the presence of massive partons. 

The relevant genuine NNLO massless and massive four-parton antenna functions re- 
quired to define our da^''^ subtraction terms for the process gg — )■ QQqq will be given in 
section 3 while their single and double unresolved limits will be given in section 4. 

Together with the mappings required to define the momenta in the reduced matrix 
elements, the subtraction terms also need a factorized form of the phase space. The 2 — )■ 
(m + 2) parton phase space needs to be factorised into an antenna phase space involving 
four-partons, the unresolved partons and the radiators, times a reduced 2 ^ m phase 
space. This factorisation is different depending on whether the hard radiators are in the 
initial or in the final state. For the cases involving only massless particles the phase space 
factorisations for final-final, initial-final and initial-initial configurations have been derived 
in [68,72,80]. The presence of massive hard radiators introduces slight differences with 
respect to the massless case in the factorisation formulae for final-final and initial-final 
configurations, but the formula in the initial-initial configuration remains unaltered as 
the radiators are massless. These factorisation formulae involving massive particles were 
derived in the final- final case in [79] and in the initial- final case in [69]. In this way, the 
integrated four-parton antennae are obtained, in all three configurations, as the integral 
over the appropriate unresolved part of the phase space of the unintegrated forms of the 
antenna functions. 



2.4 Colour-unconnected unresolved radiation subtraction: d^^j\fi^o 

In this configuration, the two unresolved partons are pair-wise colour-unconnected. The 
subtraction terms are made of two disjoint three-parton tree-level antennae of the types 
-^ijk -^ijk (^"^) "^ikj (^"^)' which have no common momenta. Together with 

those, the subtraction terms contain the colour-ordered reduced (m -|- 2)-parton amplitude 
squared |A^TO+2p (with two partons less than the original matrix element squared) and the 
jet function jin^\ The reduced matrix element squared and the jet function, both involve 
remapped momenta which arise from two independent 3—7-2 mappings appropriate to the 
given configuration. The explicit expressions for this subtraction term have been given 
in [69]. 
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This term captures the double unresolved behaviour associated to two colour-unconnected 
unresolved partons in the real matrix element squared and compensates spurious double 
unresolved features arising in subtraction terms for configurations (a). 

To obtain the integrated form of this counterterm, the (m + 2)-parton phase space 
is factorised into an m-parton phase space multiplied by two independent phase space 
factors for each of the two antennae. The integrated form is thus simply the product of 
two integrated three-parton antennae. 



3. Antenna functions 

In this section, we present the new massive four-parton antenna functions which are needed 
to define the subtraction terms for the double real radiation corrections to hadronic heavy 
quark pair production due to the process gg — )■ QQqq and presented in section 5. 

Antenna functions are the key ingredients needed to build subtraction terms in the 
antenna subtraction method and can also be used as evolution kernels in parton show- 
ers [83-86]. In general, these functions are denoted with the character X. Each antenna is 
determined by its particle content (hard radiators and unresolved particles) as well as by 
the pair of hard particles that it collapses onto in its singular limits. Antennae that collapse 
onto a quark-antiquark pair are X = A for qgq and qggq, X = B for qq'q'q and X = C for 
qqqq. Antennae that collapse onto a (anti) quark and a gluon are X = D for qgg and qggg, 
and X = E ioT qq'q' and qq'itg. Finally gluon-gluon antennae are X = F ion ggg and gggg, 
X = G for gqq and gqqg, and X = H for qqq'q' ■ All antenna functions are derived from 
physical colour-ordered matrix elements squared. Each type of antenna is obtained from a 
difi^crcnt physical process: quark-antiquark antennae are related to processes of the form 
7* — )■ qq+{paictons) [87], quark-gluon antennae, to x — ^ 5-|-(partons) [70], and gluon-gluon 
antennae are derived from H ^(partons) [71]. 

The tree-level antenna functions are obtained, both in the massless or massive case, by 
normalising the colour-ordered three- (or four)-parton tree-level squared matrix elements 
to the squared matrix element for the basic two-parton process, omitting couplings and 
colour factors: 



.0 _ c l-^i'ifei'^ 



■^ijkl - ^ijklJL |^o^|2 • i-^-^J 

S denotes the symmetry factor associated with the antenna, which accounts both for po- 
tential identical particle symmetries and for the presence of more than one antenna in the 
basic two-parton process. 

The three and four-parton initial-final and initial-initial antenna functions are, in prin- 
ciple, defined by crossing one or two massless partons from the final to the initial state 
in the corresponding final-final antennae. The three-parton initial-final and initial-initial 
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antennae are denoted respectively by X^j/^ and ^{j^ j '^l^i^^ the corresponding four-parton 
antennae will be denoted by X^jf^^ and X^f, ^. 

While at NLO only three-parton tree-level antennae are used, at NNLO also four- 
parton tree-level and three-parton one-loop antenna functions need to be considered. The 
latter, however, are not needed for the double real radiation contributions, and will not be 
discussed here. 

In contrast with the massless antenna functions, massive antennae have explicit mass 
terms and they can be of two different natures: flavour-conserving and flavour-violating. 
The latter are derived from partonic processes with a flavour-violating vertex connecting 
radiators of different flavours, one of them being massless and the other massive. Massive 
three-parton flavour violating antennae have been derived in [78], and one specific B-type 
four-parton antenna of this kind has been derived in [69]. It was needed for the construction 
of the subtraction terms for the purely fermionic processes contributing to the hadronic 
production of a tt pair at NNLO. 

For the subtraction terms that will be presented in section 5 for the partonic process 
99 QQqQj the three-parton antennae that are needed are of A, D, E and F-type in 
different configurations with massive and/or massless partons. All these antennae have 
been computed and integrated in [67,72,78,80], with the exception of one flavour-violating 
A-type antenna. This antenna function, which involves a gluon in the initial state and 
massive quark and a massless antiquark in the final state will be presented together with 
its integrated form and infrared limits in the appendix A. The unintegrated form of all the 
other required three-parton antennae will also be presented there. 

The genuine NNLO four-parton antennae that are needed for the partonic process 
under consideration are of B, E and G-type. The B and G-type antennae are known and 
will be given in the appendix B. The E type antennae are new and shall be derived below. 

For the labelling of the partons in the antenna functions used throughout this paper we 
shall use the following conventions: Massless quarks will be indexed with q while massive 
ones with Q and their mass with mg. Partons crossed to the initial state are denoted with 
a hat. The first and the last particles in the argument of a given antenna are the hard 
radiators, while the partons placed in between the radiators are the unresolved particles. 
In order to make the mass-dependence in the expressions of the antenna functions explicit, 
we define our invariants as Sij = 2pj • pj. Finally, for conciseness, the 0{e) pieces of the 
antenna functions will be omitted. 

The B-type antenna required to define our subtraction terms is massive and is only 
needed in its final-final form. It is obtained from the ratio of the matrix elements squared 
for the physical processes 7* — )• QQqq and 7* — )• QQ, and it is used to subtract the infrared 
singularities associated to the emission of an unresolved qq pair between the massive QQ 
radiator pair. This final-final massive B-type antenna is known in unintegrated and inte- 
grated form [79], and will be given in the appendix for completeness. Its infrared limits 
will be recalled in section 4. 

^It is woth noting that, in the initial-final case this crossing may not be unambiguous and may result in 
the necessity of a further splitting of the antenna as noted in [72,78]. 
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We also use a massless four-parton G-type antenna in its initial-initial form. The 
expression of this antenna is obtained by crossing two gluons in the corresponding final- 
final massless G-type antenna, which has been defined in [67] as the ratio of the processes 
H ggqq and H gg. This antenna is needed to subtract the infrared singularities 
associated to the emission of a massless qq pair radiated between the initial state gluons. 
Its unintegrated form will be given in the appendix while its infrared limits, which have not 
been documented so far, will be presented in section 4. Note as well, that the integrated 
form of this antenna has just been derived in [77]. 

The new massive four-parton E-type antenna functions are needed in their initial-final 
form. Those can be obtained by crossing the gluon in the corresponding final-final E-type 
antennae. Those are derived from the ratio of the processes x ~^ 99QQ X ~^ 99 with 
the massive gluino g playing the role of the massive (anti) quark of mass niQ. The full 
amplitude for the process x ~^ 99QQ contains leading and subleading colour pieces [71]. 
By squaring the leading colour piece, in which the qq is emitted between the gluino and 
the gluon in the colour chain, the antenna is obtained, while, squaring the subleading 
colour piece, in which the gluon is emitted between the qq pair, the antenna is obtained. 

The infrared limits of these E-type antennae derived below will be given in section 4. 
Note also that, there is on-going work [88] towards the integration of these initial- final 
massive antennae. 

To account for the infrared limits associated to the emission of an unresolved qq pair 
between a massive (anti) quark and an initial state gluon we use the following antenna 
function, 

TpOn 1 A o \ ^ / ^13*14 , Sl4 r ,2,21 

^4 3g, Ag, 2g) = -2 <^ — ^ + -^13^23 + S13 + S^l 

(Q'^ + ml) U12S234 S12S34S24 
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Sl2S34'Sl34 
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where Q"^ = -(^1+^3+^4-^2)^, rn^ 



"x — vQ^' ■sm = 513+514+534 and S234 = 534-523-524- 
In order to account for the triple colhnear hmits that involve a massless qq pair and 
an initial state gluon in those sub-leading colour amplitudes in which the gluon is placed 
between the quark and the antiquark in the gluon chain, we employ the following 
antenna 
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with Q^, rriy^, S134 and S234 given as above. Both E-type antennae are normahsed to the 
tree-level two-parton matrix element 

iM'^ixg ^ = 4(1 - e) [Q2 + ml]\ (3.5) 

4. Infrared limits 

The factorisation properties of QCD tree-level squared amplitudes have been extensively 
studied in [89—92]. While at NLO only single soft and collinear singularities may arise, at 
NNLO there are several double unresolved regions of phase space in which the amplitudes 
can develop singularities when integrated over. In the first part of this section, we shall 
list all single and double unresolved factors arising in the unresolved limits of the antennae 
needed to construct our subtraction term presented in section 5, while in the second part 
we shall list the infrared limits of all the required four-parton antennae. 

4.1 Single unresolved factors 

In their infrared limits, colour-ordered matrix elements factorise into universal unresolved 
factors and reduced matrix elements. Depending on whether the partons involved in each 
limit are massless or massive, and depending also on whether one or two partons become 
unresolved, the universal factors will be different. We shall start by presenting the single 
unresolved massless and massive factors. 

Retaining only those infrared limits which lead to poles in e when integrated over 
the phase space, the only single unresolved limits that the colour-ordered amplitudes for 
99 ~^ QQqQ have, are collinear limits involving the massless final state qq pair, and initial- 
final collinear limits involving one of the initial state gluons and the massless final state 
quark or antiquark. In the latter limit the collinear partons cluster into an initial state 
(anti) quark, while in the former, the collinear quark-antiquark pair clusters to form a gluon 
in the final state. As a consequence, the reduced colour-ordered matrix clement associated 
with this limit involves a final state gluon. Thus, although only quarks are involved in the 
final state of the process that we are presently considering, unresolved gluon limits have to 
be considered as well. 

When two massless partons become collinear, a colour-ordered sub-amplitude factorises 
into a reduced matrix element times a specific Altarelli-Parisi splitting function correspond- 
ing to the particular parton-parton splitting. These functions depend on z, the momentum 
fraction carried by the unresolved parton. Depending whether the unresolved parton is 
collinear to an initial or to a final state parton, the definition of z is different. For two final 
state particles i and j of momenta Pi and pj becoming collinear, we have, in the limit, 

Pi zpij, Pj ^ (1 - z)pij, Sik zsijk, Sjk ^ (1 - z)sijk , (4.1) 

whereas for a final state particle j of momentum pj becoming collinear with an initial state 
parton i of momentum pi we have 

Pj zpi, pij ^ (1 - z)pi, Sik -f^, Sjk . (4.2) 
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The splitting functions denoted by Pij^(^ij-^{z) corresponding to the colhnear Umit of 
two final state partons i and j, given in the conventional dimensional regular isation scheme 
with all particles treated in d = 4 — 2e dimensions, are [93]: 

W^)= ^^'^-f (4.3) 

F«.„(.) = ''^';:f-- (4.4) 



F»^g(2) = 2 



Z I — Z , 

+ + z{l - z) 



(4.5) 



When the collinearity arises between an initial (i) and a final state parton (j), the 
splitting functions denoted by Pjj.(_(jj)(z) are given by [93]: 

, , l + z'^-eil-zf 1 1 , , , , 

1 + fl — zl^ — fZ^ 1 
P,9^Qi^) = _ = Y^/,9^Qi^) (4.7) 

+ fl — — f 1 — f 
= = ^^P,,^g{z) (4.8) 

2(1- z + z^)'^ 1 
^..-g(.) = = y^/99->g{z). (4.9) 

The additional factors (1 — e) and 1/(1 — e) account for the different number of polarizations 
of quark and gluons in the cases in which the particle entering the hard processes changes 
its type. 

In all the splitting functions defined above, the label q can stand for a massless quark or 
an antiquark since charge conjugation implies that Pqg^Q = Pqg^Q ^iid Pqg-i^Q = Pqg^Q ■ 
The labels Q and G denote the parent parton of the two collinear partons, which is massless. 



When a gluon j emitted between the hard radiators i and A; becomes soft, the eikonal 
factor that factorises off the colour-ordered squared matrix element is 

where rrii and m^ arc the masses of the radiators. This limit is obtained by letting pj — t- \pj 
with A — )• or, alternatively, by keeping only terms with any product of the inverse of 
the two invariants Sij and Sjk in the colour-ordered matrix element squared. By letting 
mi,mk the usual expression for the massless soft eikonal factor is obtained. 



4.2 Double unresolved factors 

In general, when two particles are unresolved in a tree-level process, a variety of different 
configurations can arise: 
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(1) two soft particles, 

(2) two pairs of collinear particles, 

(3) three collinear particles, 

(4) one soft and two collinear. 

As we saw in section 2, the resulting double unresolved configurations of the colour- 
ordered matrix element squared need to be separated into three categories depending on 
the colour connections of the unresolved partons and the hard radiators associated with 
them. In the following, we shall describe only the double unresolved factors encountered 
in the context of this paper, which correspond to the configurations (1), (2) and (3) listed 
above, but not to (4), since the process gg — )■ QQqq does not have any soft-and-collinear 
limits. 

4.2.1 Colour-unconnected pairs of unresolved particles 

For unresolved particles that are disjoint in the colour chain, which arise in the item (d) 
as presented in section 2, the colour-ordered matrix elements squared factorise into the 
product of two disjoint single unresolved factors nuiltiplicd by a reduced matrix element 
with two partons less than the original colour-ordered matrix element squared. 

For example, in the case of two pairs of colour-connected particles (a, b) and (c, d) 
which are all located in the final state, the colour-ordered matrix elements given by 

(. . . , a, 6, . . . , c, d, . . .)p undergo the following factorisation in the double collinear limit: 

|A^n(- . . , a, 6, . . . , C, d, . . .)P ^ Pab^p{zi, Sab) Pcd^Q{z2, Scd)|>ln-2(- • • , ■ ■ ■ , Q, • • 

(4.11) 

where partons a and b form P, while c and d cluster to form Q, so that P and Q are 
themselves colour-unconnected. The collinear splitting functions appearing in this equation 
are the (colourless) Altarelli-Parisi splitting functions given by eqs. (4.3-4.9). 

In the process that we are considering in this paper, the only double unresolved colour- 
unconnected limits that can occur are double collinear (anti) quark-gluon limits, with each 
single collinear pair given by a final state (anti) quark and an initial state gluon. The 
unresolved factor associated to this limit is a product of two splitting functions of the 
type given in eqs. (4.3-4.9). Note that none of the four-parton antennae required for our 
subtraction term captures these double collinear singularities. These are accounted for by 
subtraction terms involving the product of two three-parton antennae. 

4.2.2 Colour-connected pairs of unresolved particles 

The colour-connected double unresolved limits the partonic process gg QQqq develops 
are: 

A) Triple collinear limits involving an initial state gluon and a massless final state quark- 
antiquark pair. 
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B) Double soft limits of a massless final state qq pair emitted between massive or massless 
radiators. 

As mentioned above, these colour-connected double unresolved limits are captured by four- 
parton antenna functions and arise in subtraction terms dcr^^^Q as presented in section 
2. In the following, we shall give the massive and massless double unresolved factors asso- 
ciated with these two types of limits. 

A) Massless triple coUinear factors 

In those regions of phase space where three colour-connected massless partons (a, 6, c) be- 
come collinear, a generic colour-ordered amplitude squared denoted by |A^Ji(. . . , a, 6, c, . . .)p 
factor ises as: 

(. . . , a, 6, c, . . .)!' ^ Pabc^plMl.^i. . . , P, . . .)!'• (4.12) 

where the three colour-connected final state particles (a, b, c) cluster to form a single parent 
particle P. The triple collinear splitting function for partons a, b and c clustering to form 
the parent parton P is generically denoted by, 

Pabc^p{w, Sabc), (4.13) 

where w, x and y are the momentum fractions of the clustered partons, 

Pa = wpp, pb = xpp, pc = ypp, with w + x + y = l. (4.14) 

In addition to its dependence on the momentum fractions carried by the clustering partons, 
the splitting function also depends on the invariant masses of parton-parton pairs and the 
invariant mass of the whole cluster. The explicit forms of the triple collinear splitting 
functions Pabc-^p are obtained by retaining terms in the colour-ordered matrix element 
squared that possess two of the 'small' denominators Sab, Sac, Sbc and Sabc- 

The triple collinear limits which have to be considered here are those involving the 
massless final state quark-antiquark pair and one of the initial state gluons. The splitting 
functions for these types of (initial-final) triple collinear limits can be obtained from the 
analoguous limit where the three collinear particles are in the final state [89,92]. 

The clustering of a gluon with a quark-antiquark pair into a parent gluon has two 
distinct functions depending the colour connection of the collinear particles. In leading 
colour contributions where the gluon is emitted "outside" the quark-antiquark pair, the 
following non-abelian splitting function is obtained 

Pgqq^oif^i Hi ^gqi ^qq-i ^gqq) ~ 

1 f A^gq I /'I ^\^qq i /o ^\\ "^(."^^gqq (1 ^)^9'?) 




^gq J ^qq^gqqi^ ^ 

1 ( i'^-y) _y_2w-^-'^^^^~y^^y~^ 



(1 — e)w{l — w) 

^) 2x^1 - yw - {1 - y)il - w)) 
w{l — w) (1 — e)w{l — w) 
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_^ 1 / + Axw) ^ 2x {w{x - y) - y{l + w)) \ 
^qq^gqq V ^^(1 ~ (1 " - w) ) ' 

In those sub-leading colour pieces where the gluon is emitted "between" the quark-antiquark 
pair in the colour chain, one obtains a QED-like splitting function 

Pqgq^dw, X, y, Sqg, Sgq, Sqq, Sqgq) = 

"qgq \ "qg J "gq^qg \ it 

1 + 2x + e - \ _ Y + {sqg o Sgq,w^ y). (4.16) 

^qg^qgq \ \^ ^) J 

The triple collinear splitting functions given in eqs.(4.15) and (4.16) correspond to 

configurations in which all three collinear particles are outgoing. As mentioned above, in 

the present calculation we do not deal with these limits (indeed, there are not enough 

massless particles in the final state to have final-final triple collinear limits). Instead, we 

deal with triple collinear limits where one of the partons (a gluon) is in the initial state. 

In this case, initial-final triple collinear splitting functions arise. They are related to their 

final-final counterparts as follows [92] 

Pgqq'k-oi^^i -I Zl, Sgq, Sqq, Sgqq) = 

Pgqq-^G0-/Z3, -Z2/Z3, -Z1/Z3, -Sgq, Sqq, Sqg - Sgq - Sqg) 
Pqgq-ir-G{zi, Zj„Z2, Sqg, Sgq, Sqq, Sqgq) = 

Pqgq-^oi, Z\/ Z"}, 1/ ^3, Z2/Z3, Sqg, Sgq, Sqq, Sqq Sgq Sqg), 

(4.17) 

where zi and Z2 are the momentum fractions final state quark and antiquark respectively, 
and zs = 1 — zi — Z2- 



B) Massive and massless double soft factors 

When a massless quark-antiquark pair becomes soft between two hard radiators, sub- 
amplitudes squared factorise into a double soft factor and a reduced matrix element squared 
with the quark-antiquark pair removed from it. As it was the case for the single soft gluon 
cikonal factor, this double soft factor, needed for a soft qq pair, also depends on the masses 
of the hard radiators if these are massive. 

When the quark-antiquark pair (c, d) becomes soft between the hard radiators (a, b) of 
masses ma and rrib, the massive double soft factor is given by, 

SacSad _|_ ^bc^bd 
(Sac + Sad)'^ (Sbc + Sbd)"^. 

^< (4.18) 

ScdiSac + Sad)'^ Scd(s6c + Sbd)^ 

This factor is obtained by setting pc — )■ Xpc, Pd ^Pd with A — >■ in the matrix elements 
squared and was derived in [79] as the soft qq limit of the massive final-final antenna 
B2iQ,q,q,Q) and in [69] using current algebra. Note also that, the corresponding massless 
factor [67, 89] can be obtained from this massive one by setting the masses to zero. 



Sacdb{ma,mb) 



'^{SabScd SacSbd S^cSad) ^ 2 
sLiSac + Sad) (Sbc + Sbd) 
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4.3 Single and double unresolved limits of four-parton antennae 

In their infrared limits, the four-parton antennae of B, E and G types required to compute 
the double real contributions to heavy quark pair production in hadronic collisions due to 
the partonic process gg QQqq, and defined in section 3, yield the universal single and 
double unresolved factors defined above. Additionally, some of these four-parton antennae 
yield angular correlation terms in single coUinear limits. These angular terms arise when a 
gluon splits into a quark-antiquark pair or into two gluons, and the way in which they are 
dealt with in the antenna subtraction method has been explained in detail in [68,69,94,95] 
and will not be recalled here. 

In addition to the conventions already defined in section 3 for the labelling of the par- 
tons present in antenna function, we shall denoted with {ij)a the momentum of the parent 
partons in collinear limits. Thus, partons labelled with {ij)a will have its momentum 
Pi + pj in final-final collinear limits, whereas its momentum will be given by pi — pj in the 
initial-final case. 

Massive final-final B-type antennae 

The massive final-final B-type antenna given in eq.(B.l) is only singular when the massless 
qq pair is soft or collinear. The behaviour of this antenna in each of these limits is given 

by 

B0(1q,4,-,3„2q) 5i342(mQ,mQ) (4.19) 

S2(1q,4,-,3„2q) —Pgg^G{z)Al{lQ,i34)g,2Q) + s,ng., (4.20) 

where (ang.) indicates the presence of angular correlation terms. 

Massless initial-initial G-type antennae 

The massless initial-initial G-type antenna given in the appendix in eq.(B.3) has the fol- 
lowing infrared limits 

^^(i^, 3„ 4g, 2g) 5i342(0, 0) (4.21) 

Gl{lg,3q,Aq,2g) > Pgqq<r-G{z3,Z2,Zl,Si3, 834,8134) (4.22) 

r, - 2g\\3q\\4g 
G^{lg,3q,'ig,2g) > -Pg^^^G ( ^3 , ^2 , ^1 , S 24 , S34 , S 234 ) (4.23) 

Gliig, 3q, Aq, 2g) "'-^ —Pqq^G{z)Fl{ig, (34),, 2,) + aug. (4.24) 
S34 

Gl{ig,3qA-qX) —Pqq^G{,z)Gl{2gAq,^^S)q) (4.25) 
Sl3 

Gl{lg,3qAq,^g) ^"-^ ^Pqq^G{z)Gl{ig,3q,{2A) ). (4.26) 
S24 ^ 

Massive intial-final E and E -type antennae 

The initial-final E-type antenna given in eq.(3.3) has a soft qq limit, a triple collinear limit. 
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a final-final and an initial-final single collinear limits. The behaviour of this antenna in 



each of these limits is 

E2{lQ,3g, 2g) 5i342(mQ, 0) (4.27) 

r. 2gl|3ql|4q 

Ed^Q,^qAq,^g) > Pgqq^G{z3, Z2, Zl, 824, 834, S234) (4.28) 

E^^ilQ,3g,4g,2g) — Pg^^G {z) DI{1q , (34) g , 2 g) + aug. (4.29) 
S34 

El{lQ,3g,4qA) —Pqq^Giz)El{lQ,3g,{24)). (4.30) 



The initial-final massive antenna D^{Q, g, g) is obtained by crossing to the initial state one 
of the gluons in its final-final counterpart. As it can be seen in eq.(4.29), it is the antenna 
onto which E^{Q, q, q, g) collapses to in its single q\ \q collinear limit. It will be given in the 
appendix A. 

The sub-leading colour antenna E^ has a QED-like triple collinear and two initial-final 
single collinear limits: 

~„ ^ 2g||3,||4j _ 

E4{lQ,3g,Ag,2g) > Pggg^G (^1 , ^3 , ^2 , S 23 , S 24 , S34 , ^234 ) , (4.31) 

^40(1q,3„4,-,2,) —Pgg^G{z)ElilQ,3g,(24) ) (4.32) 

S24 ^ 

^0(1q,3„4,-,2,) —Pgg^G{z)El{lQ,Ag, (23)^). (4.33) 

S23 

5. Top quark pair production at the LHC 

In this section we shall present the double real emission contributions to ti production 
at the LHC due to the process gg — ^ QQqq. Together with these, we shall give their 
corresponding antenna subtraction terms, which capture all single and double unresolved 
limits of the leading and subleading colour pieces of the real radiation matrix elements 
squared. 

5.1 Conventions 

To facilitate the reading of our expressions, we shall closely follow the notation in [69, 78] 
for matrix elements and subtractions terms. The main points of our conventions are the 
following: The matrix elements denoted with M represent colour-ordered sub-amplitudes 
in which the coupling constants and colour factors are omitted. Furthermore, to explicitly 
visualise the colour connection between particles in these colour-ordered amplitudes, a 
double semicolon is used in the labeling of the partons present in a given matrix element. 
This double semicolon is used for separating chains of colour-connected partons. Partons 
within a pair of double semicolons belong to a same colour chain, and adjacent partons 
within a colour chain are colour-connected. An antiquark (or an initial state quark) at the 
end of a colour chain and a like flavour quark (or initial state antiquark) at the beginning of 
a different colour chain are also colour-connected since the two chains merge in the collinear 
limit where the qq clusters into a gluon. We also denote gluons which are photon-like and 
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only couple to quark lines with the index 7 instead of g. In sub-amplitudes where all 
gluons are photon-like no semicolons are used, since the concept of colour connection is 
not meaningful. A hat over the label of a certain parton indicates that it is an initial state 
particle (for example, Iq is an initial state quark with momentum pi). 

Concerning the notation in the subtraction terms, the conventions for the reduced 
matrix elements are the same as those for the real radiation matrix elements discussed 
above. The remapped final-state momenta are denoted with tildes and the remapped mo- 
menta of initial state hard radiators are denoted by a bar and a hat, as used in other 
papers [68,69,75]. In the four-parton antenna functions the hard radiators are "on the 
edges" and the uresolved particles are "in the middle" . 

5.2 Double real radiation contributions 

For two incoming hadrons, Hi,H2, the hadronic heavy quark pair production cross section 
may be written as 

do- = ^y^ fa/li^ljfJ-p) fb/2i^2,fJ'F)daab{^lHi,^2H2,HF,fJ'R) ■ (5.1) 

a,b 

^1 and ^2 are the momentum fractions of the partons of species a and b in both incom- 
ing hadrons, fi being the corresponding parton distribution functions, daab denotes the 
parton-level scattering cross section for incoming partons a and b which depends on the 
the renormalisation and factorisation scales denoted by and iip respectively. The par- 
tonic cross section d&ab has a perturbative expansion in the strong coupling ag which itself 
depends on the renormalisation scale fiR. 

Following the general factorisation formula given in eq.(5.1) (if we omit the renormalisa- 
tion and factorisation scale dependences) the contribution of the partonic process gg QQ 
to the leading order cross section for heavy QQ production in hadronic collisions is 

The leading order partonic differential cross section, written in terms of colour-ordered 
matrix elements, is given by 



d&gg_^QQ = ML0d^2 {Pl ,P2;P3,P4) 

1 



7V4|>tO(lQ, 3„ 4„ 2q)|2 + |>tO(lQ, 4„ 3„ 2q)|2) 

(pi,P2). (5.3) 



In this equation, the normalisation factor Mlo reads 

where s is the hadronic center of mass energy, {N^ — 1) comes from the colour sum, the 
factor 1/4[N^ — 1)^ accounts for the averaging over the spin and colour of the incoming 
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two gluons, and (l/2s) is the hadron-hadron flux factor. The couphng is deflncd as usuah 
as = g'^/in, C(e) = (47r)^e"^'^/87r^, and (^(e) = {AnYe^^'^. This way of expressing the 
couphng factors, with each power of ag accompanied by a power of C'(e), keeps the couphng 
dimensionless in dimensional regularisation. The phase space for the production of a QQ 
pair of momenta pi and P2 is given by d^2{pi,P2',P3,P4)- In it, ps and P4 are the momenta 

(2) 

of the initial state gluons. The jet function denoted by J2 ensures that the top and the 
antitop are in two separate jets. 

The colour-ordered amplitudes Ai^ are related to the full amplitude M4 (Ig, 2q, 3g, 4^) 
through the colour decomposition 



M|(1q,2q,3„4,) = g\V2f\^{T-^T-X,^Ml{lQ,3g,lg,2Q) 

+ {T^'T^%,,Ml{lQ,lg,3g,2Q)y (5.5) 

and the colour-ordered amplitude with photon-like gluons is given by 

MlilQ, 3^, 4^, 2q) = Ml{lQ,3g, Ig, 2q) + Ml{lQAg, 3„ 2q). (5.6) 

The contribution of the gg — > QQqq partonic channel to the double real radiation cross 
section for heavy quark pair production in hadronic collisions is 

with the partonic cross section given by 

d^fg^QQqq = ^NNLO d^4 (pi , P2 , P3 , P4 ; P5 , P6 ) | ( 1q , 2q , 3^ , 4g , 5g , 6g ) | ^ ^ (pi , P2 , P3 , P4 ) ■ 

(5.8) 

The full matrix element squared |Mg (. . .)p is summed over spin and colour. The normal- 
isation factor A/jy^^Q accounts for the spin and colour averaging, the flux factor and the 
sum over the possible flavours of the massless qq pair. The jet function j'^'* corresponds 
to the selection criteria of a 2-jet event: out of four partons, from which two are a QQ 
pair, an event with two jets is built. Each of these two jets has the heavy quark Q or the 
heavy antiquark Q in it. The additional partons present in a jet are either theoretically 
unresolved (soft or collinear) or not "seen" by the experimental resolution criteria. 

The full matrix element for the process gg QQqq has the following colour decompo- 
sition [78] 

<(lQ,2g,3„4,-,5„6,) = 

{T^iT^i),^,^5i,,i^Ml{lQ:ig3g. %; ; 3„ 2q) 



9\v2f y: 

{i,j)eP(5,6) 



+ (r«0ni4(r"^)i3i2K(lQ>^9>4,-; ; 3„i„ 2q) 
+<5i,,i,(T"*T"0i3i2-^6(lQ, %; ; 3g, ig,j9^ 2q) 



1 



-_(^r''iT-^)i^iJi,,i,MtilQ,igJg,2Q-,3g,Ag) 
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(5.9) 



Squaring eq.(5.9), summing and averaging over spin, colour, and quark flavour, and plug- 
ging our result in eq.(5.8) allows us to write the partonic double real radiation cross section 
in the following form 



a,\2 C(e)2 , 

1 ^T^d$4(Pl,P2,P3,P4;P5,P6) 



^ ' {jj)eP(5,6) '- 

' (|A^[!(1q, 45; ; 3„ 2q)|2 + |A^0(1q, i„ 4,; ; ^^jg, 2q)|2 



+|7W°(lQ,4,-;;3„i^,jg,2Q)|2) 

+\Ml{lQ,ig3g. 2q; ; 3„ 4,-)|2 + \mI{\q, ig, 2q- ■ 3„ig, 4,-)|2 
+|7W°(1q, 2q; ; 3„ ig,jg, 4g-)|2 - |A^0(1q, J^, j^, 4g-; ; 3,, 2^)1^ 
-|>fO(lQ, ig, 4,-; ; 3„ig, 2q)|2 - \MI{Iq, 4,-; ; 3<„ 2q)\^ 
+2Re(A^0(lQ, 4,-; ; 3„ 2q)MI{Iq, 4,-; ; 3^, Jg,]^, 2^)^^) 
+2-Re{Ml{lQ,ig,jg, ; 3,, 2q)MI{Iq, 4,-; ; 3„ jg, i^, 2q)^) 
+Re{Ml{lQ,ig, 4q-; ; 3q, jg, 2Q)Ml{lQ,jg, 4,; ; 3^, Zg, 2^)^) 
-Re(7W|](lQ, Zg, jg, 4g-; ; 3g, 2q)>1°(1q, jg, Jg, 4g-; ; 3g, 2Q)t) 

-Re(A^°(lQ, 4g-; ; 3g, igX 2q)MI{Iq, 4g-; ; 3g, jg, ig, 2Q)t) 

-2Re(7W[|(lQ,ig,ig,4,-;;3g,2Q)A^0(lQ,ig,ig,2Q;;3g,4g-)t) 



-2Re(7W[|(lQ,ig,4g;;3g,jg 
-2Re(7W^(lQ,4g-;;3g,ig,jg 
-2Re(7W^(lQ,ig,jg,45;;3, 
-2'Rc{Ml{lQ,igAqV,^q,jg 
-2Re{Ml{lQ,igAqV,^q,3g 
-2Re(A^[|(lQ,4g;;35,^g,jg 
-2Re(7W0(lQ,ig,i3,4,^;;3, 
-2Re(7W[|(lQ,4,-;;3y,ig,jg 
-2Re(7Wg(lQ,ig,4^;;3g,jg 
1 /l..,o 



2Q)A^0(lQ,ig,ig,2Q;;3„4,-)t) 
2Q)A^0(lQ,?g,ig,2Q;;3„4,-)t) 
2Q)A^S](lQ,ig,2Q;;3„jg,4,-)t) 
2Q)A^[](lQ,ig,2Q;;3„jg,4,-)t) 
2Q)A^0(lQ,jg,2Q;;3„ig,4,-)t) 
2Q)A^^(lQ,J,,2Q;;3„ig,45)t) 
2Q)A^0(lQ,,2Q;;3„ig,ig,4,-)t) 



2q)A^0(1q,,2q 
2q)K(1q,,2q 

2 I A/fO 



:3q.%.J,g,45)^) 

j35,Jg,ig,4^)^) 



+ iV2 (^2 '-^^(Iq, 2q, 3„ 4,-, i^,i^)|^ - |>t^(lQ, zg,ig, 2q; ; 3g, 4g-)|^ 

-|A^0(1q, ig, 2q; ; 3g,ig, 45)|2 - \MI{Iq,2q- ; 3„ ig,Jg, 4,-)|2 
+Re(7W[!(lQ, ig, 2q; ; 3,,ig, 45)7W[!(1q, jg, 2^; ; 3^, ig, 4g-)t) 
-Re(7W[!(lQ, ig,ig, 2q; ; 3^, 45)7W[!(1q, jg, ig, 2q; ; 3,, 4g-)t) 
-Re(7W[!(lQ, 2q; ; 3^, ig, jg, 4g-)7W[!(lQ, 2q- ; 3^, jg, ig, 4g-)t) 
+2Re(A^0(lQ, ig, jg, 2q; ; 3,, 45-)A^0(1q, 2^; ; 3,, ig, jg, 4,-)t) 
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4^^(Pl,P2,P3,P4)- 
(5.10) 



As it can be seen in eq.(5.10), this cross section does not only contain colour-ordered 
matrix elements squared, but also interference terms of two different colour-ordered matrix 
elements. 

5.3 Subtraction terms 

In this subsection we present the subtraction terms which capture all unresolved behaviour 
present in the real contributions given above in eq.(5.10) and related to the partonic process 
99 QQqq- As explained in section 2, the full subtraction term receives in this case three 
different contributions 

da^ ^r^r, - = da^'"' _ + da'^''' _ + da'^''^ _. (5.11) 
99-^QQqq gg^QQqq gg^QQqq gg^QQqq ^ ' 



The da*^'" terms subtract the single unresolved limits of the six-parton real radiation matrix 
elements. They are constructed as products of three-parton trcc-lcvcl antenna functions 
and fivc-parton reduced matrix elements with remapped momenta. The do"^''' pieces arc 
genuine NNLO subtraction terms which account for those double unresolved limits of the 
real radiation matrix element that involve a pair of colour-connected partons; in this case 
these can be double soft and triple collinear limits. Finally, the da^''^ terms, constructed 
as a product of two three-parton tree-level antenna functions and a four-parton reduced 
matrix element, capture the double collinear behaviour of the real matrix clement squared 
and compensate for the over-subtraction of colour-unconnected double unresolved limits 
introduced in do"^'". In the following, we shall outline in some detail the derivation of 
do"^'^ for which we follow a different approach than the one usually employed in the con- 
struction of antenna subtraction terms. Our starting point is the factorisation of colour 
sub-amplitudes in their soft qq limits. 

5.3.1 Double soft behaviour of amplitudes 

As it was explained in [69] when a soft quark-antiquark pair is emitted between partons a 
and h in the colour chain, the colour-ordered amplitude denoted by (..., a, Cq\ ; d^, 6, ...) 
factorises as 

■A^n(-,a,Cq; ...) ^ \usJ<Pd)liiVsXVc)\ {JaiPcPd) - JbiPcPd)) Ml_2{:;a,b, ...), 

(5.12) 

where the soft currents are given by 

1 

This factorisation can be applied to the first three lines of eq.(5.9). Squaring eq.(5.12) and 
summing over the spins of the soft quark and antiquark gives 

\Ml{...,a,Cq;;dq,b,...)\^ Sacdbima,mb)\Ml_2{...,a,b, ...)\^ (5.14) 



JiiPJ^Pk) = „ \ „ ^ - (5.13) 
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with the double soft factor Sacdbi^ai'nih) given by 
Sacdb{ma,mb) = tr(i)d7^^c7i/) {JaiPcPd) - JbiPcPd)) (JaiPcPd) - JbiPcPd))^ ■ (5-15) 

By expUcitly evaluating the trace in eq.(5.15) and using the definition of the currents in 
eq.(5.13) the expression given in eq.(4.18) for the double soft factor is obtained. 

When a qq pair becomes soft in a sub-amplitude whose colour factor contains Si^^i^, like 
MQ{lQ,ig, jg,2Q; ;3q,Aq) in the fourth line of eq.(5.9), the factorisation given in eq.(5.12) 
does not hold. This is due to the fact that the presence of the factor indicates that the 
gluon propagator from which the soft pair splits is photon-like. In these cases, factorisation 
at the amplitude level reads 

Ml{...,a;;dq,Cq)'"'^ [usd{Pd)l^lVpSPc)]i ^ Jf{Pc,Pd)- ^ -/f (PcPd) J>tn-2(->a) 

(5.16) 

where {q} is the set of all final state quarks and initial state antiquarks in the partonic 
process, and {q\ is the set of all final state antiquarks and initial state quarks. In the 
process gg — > QQqq, however, since there is only one quark-antiquark pair in addition to 
the soft one, each sum in eq.(5.16) contains only one term. Also, the sub-amplitudes in 
the last two lines of eq.(5.9) are finite in the soft qq limit as the quark and the antiquark 
are not colour connected. 

With these considerations, we can obtain the double soft behaviour of the full amplitude 
for the process gg — > QQqq by readily applying eqs.(5.12,5.16) to eq.(5.9). It reads, 

M0(lQ,2Q,3„4,-,53,6g) g\^f\usMli^'^s,{PA)\ 

{i,j)eP{5,6) ^ 

HT''')nu{T''')isi.{4iP4,P3) - J^{P4,P3)) 

-^^{T'''T-')niAs,u{jliP4,Ps) - J^{P4,P3)) M'i{lQ,ig,jg,2Q). 

(5.17) 

We can now square eq.(5.17) and evaluate the Dirac traces as well as the colour traces to 
obtain the double soft limit of the full matrix element squared in terms of reduced colour- 
ordered matrix elements squared (with two partons less than the original amplitude) and 
combinations of double soft factors: 

\M^{1q, 2q, Sq, 4q, 5g, 6g)f g^N^ - 1) 



143i("T'Q,0) +Sj432 (0,mQ)+5,43,(0,0))|A^0(lQ,z„j„2, 
-Su32{mQ,mQ)(\Ml{lQ,ig,jg,2Q)\^ - {l/2)\MlilQ,i^,j^,2Q)f 
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■ (Su3iimQ, 0) + iSj432(0, mQ^j \MlilQ,ij,jj, 2^ 



■^^Su32imQ,mQ)\M^{lQ,i^,j^,2Q)\ 



(5.18) 



In eq.(5.18), we see that in the soft factors the hard radiators are two particles taken from 
the hst {1q, 2Q,ig,jg}, and depending on which of them are involved (massive or massless) 
the soft factors will contain two, one, or no mass terms. 

5.4 The construction of dfj^ '' „^ _ 

gg^QQqq 

We decompose the da^'^ subtraction term into a part which can be directly related to 
the double soft limit of the amplitude squared which we denote do-'^'''^^^ and another that 
accounts for all other double unresolved behaviour of the real matrix element for the process 
99 QQqQ which is not accounted for in da^'^^^^ . We shall denote this latter part of the 
subtraction term as da'^'^'^^^ such that, 

d(T - = da _ + dcj _. (5.19) 

gg->-QQqq gg-^QQqq gg->-QQqq ^ ' 

From cqs. (4. 19,4.21,4.27) we know that 5i432(mQ,mQ) is the double soft factor present 
in the double soft limit of B2{lQ,4:q,3q,2Q), Si4Si{mQ,0) appears in the double soft limit 
of E2{lQ,3q,4:g,iy), S243i{mQ,0) appears in the double soft limit of £'4(2^, 4g, 3g, ip), and 
finally 5i43j(0, 0) is the double soft limit of Gl{ig, 3q,4q,jg). 

The subtraction term da^'^^^^ is therefore obtained by replacing in eq.(5.18) the double 
soft factors with the corresponding four-parton antennae, remapping the momenta in the 
reduced matrix elements accordingly, and removing all spurious single collinear limits of 
the four-parton antennae with terms of the form ■ X^. It reads, 

'^^g^Qqq = ■^^oNf (g:) ^^dMPUP2,P3,P4;P5,P6) Yl 



+ 



(^Ei{lQ,3qAqJg)-G'i{ig,3q,Aq)Dl{lQ,{34)g;ig) 

-^°(lQ,?,,4,-)E30((M)Q,3,,l,))|>t°((m)Q,t,,i,,2Q)|2jf (^4,P2) 

(^El{2QAq,3q^g)-Gl{ig,3q,Aq)Dl{2Q,{34)g;ig) 

-Al{2Q:ig,3q)Eli{23)Q,4q,tq))\Mli^^ 

+ (Gi{ig,3q,Aq,Jg) - {l/2)G'i{ig,3qAq)F§{ig,mg,jg) 

-il/2)Gl{jg,^Aq)F^(ig, mgJg) ' a°(3„ 4,-)G°(j„ (3~4),-,l,-) 
-a0(4,-,jg,3,)G0(z\,(3l)g,jg))|A^0(iQ,i,,j5,2Q)|Vf (pl,P2) 

- (£4°(1q, 3„ 4q, ig) - Gliig, 3q, 4q)Dl{lQ, (3~4)„l,) 

-Al{lQ,ig,4q)Eii{U)Q,3qX))\Mli{mQ,t^J^,^^^^^ 
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-(^4(2Q,%,3„i,) - GO(i„3„%)i?3°(2o, (34)5,^,) 

- (bO(1q, 4,-, 3„ 2q) - {l/2)El{lQ,3q, 4q)Al{{lA)Q, (34),, 2q) 

-(1/2)^3 (2q,3„4,-)^0(1q, (34),, (21)q)) x (|>10((134)q, i„ j„ (234)q|^ 
-(l/2)|M5((134)Q,i^,i^,(234)^|2)jf(pir4,p^) 

(52(1q,4,-,3„2q) - (l/2)£;30(lQ,35,45-)^0((r4)Q,(3'4)„2Q) 

-(1/2)^30(2^,3,, 4,-)A0(1q, (3~4)„ (2~4)q)) 



iV2 



:|A^^((134)Q,i,„jV,(234)Q|2jf)(pir4,P^4 



(5.20) 



As explained in section 4, the initial-final massive antenna appearing in this sub- 
traction term, is the "full" antenna. It is this antenna onto which the antenna function 
collapses to in its single qq collinear limit. 

In order to obtain the full da^'^ subtraction term, we have to supplement eq.(5.20) with 
additional terms that ensure the correct subtraction of all other colour-connected double 
unresolved limit of the real radiation matrix clement squared. The only other unresolved 
limits of this type that the process gg — )• QQqq has are triple collinear limits ig||3q||4q. In 
these limits the full matrix element squared factorises as [89-92] 



|mO(1q,2q,3„4,-,5„6,)| 



9 'sl|3gl|4<, 



{Pgiqiqs-h-G + Pq4q39i^G) 



qsQiqA-^G 



|M°(lQ,2Q,(34z)^,j,)|2. (5.21) 



Taking now the triple collinear limit of da^'^^^^ we obtain, 



^""qg-^QQm ~^ ^'''^^^{2^) ^(^d$4(pi,P2,;>3,P4;P5,P6) 



PqAqsgi^G 

(5.22) 

where we have used eqs. (5.5,5.6) to relate the full four-parton reduced matrix element 
squared denoted by IM4 (Ig, 2q, (34z)^, j,)p to the corresponding colour-ordered ampli- 
tudes. The triple collinear splitting functions Pijk^-G have been defined in section 4. For 
conciseness, in the above equations the arguments of these splitting functions have been 
omitted. 

We see that the subtraction term da^''^^^\ originally derived in order to capture the 
soft qq limit, also subtracts the non-abelian piece of the triple collinear limits i,||3g||4g. To 
account for the QED-like triple collinear limits, proportional to 1/Nc in eq.(5.21), we use 



-29- 



initial-final massive four-parton antennae. The following subtraction term is obtained: 

-^(^4°(lQ,3g,%,g-A0(lQ,i,,4,-)S°((r4)Q,3,,l,) 

-A'i{lQ;ig,3g)El{{r3)Q,^,tg)) X (|7W2((134)Q,i5,j^,2Q|2 
+|M^((134)Q,i^,t5,2Q|2)jf (pir4,i?2) 
-^(^0(2Q,4,-,3,,g-A0(2Q,i,,3,)^3°((23)Q,4,-,l5) 

-Al{2Q:ig,^)El{{M)Q,3j,)) X (|>lO(lQ,i„j„(234)Q|2 
i„ t„ (234)q (pi , p^4) 
^(^40(lQ,3,,4,-,z,)-AO(lQ,i„4,-)^?30((r4)Q,3,,?,) 

-Al{lQ;ig,3q)Ei{{13)Q,4g,ig))\Ml{{lM)QXh^ 
+ ^(^4(2Q,4,-,3„g-A0(2Q,i„3,)^;°((2'3)Q,4,-,i,-) 

-^°(2Q,^^,4,-)£;3^(2l)Q,3,,l,))|>t^(lQ,l^,j^,(234)Q|Vf^ 



+ 

1 



(5.23) 

Since the £"4 antenna docs not possess any soft qq limit, the subtraction term given in 
eq.(5.23) does not have any soft qq limits either and only subtracts the abelian triple 
collinear limit, as required. The sum da^'^^^^^ + da^'''^^^ will therefore correctly subtract 
both double soft and triple collinear limits of the double real radiation matrix element 
squared, without introducing any spurious single unresolved singularities. 

5.5 Construction of the do"^'" „^ _ and da^''^ „^ _ subtraction terms 

gg^QQqq gg-^QQqq 

The do"^'" subtraction terms are NLO like: They are constructed as products of three- 
parton tree-level antennae and five-parton reduced matrix elements, and they subtract the 
single unresolved limits of the real radiation matrix elements. Thus, for the present calcu- 
lation, they could a priori be taken over from [78], where we derived the NLO subtraction 
terms for tt + jet production. 

However, at NNLO, the five-parton reduced matrix elements with remapped momenta 
present in these NLO-like subtraction terms develop further single unresolved limits. Those 
singularities are unphysical, since they do not correspond to any unresolved behaviour of 
the real radiation matrix element squared and must be cancelled by the X^-X^ pieces of the 
b-type subtraction terms. In dcr"^''', these X^-X^ terms are dictated by the requirement that 
all single unresolved limits of the four-parton antennae should be removed. As the choice of 
these four-parton antennae is fixed by the double unresolved behaviour of the real radiation 
matrix element, we do not have the freedom to choose the three-parton antenna functions 
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that make do"^'" simplest, as we did in [78]. Instead, these three-parton antennae have to 
be carefully chosen in such a way that the cancelation of the aforementioned unphysical 
singularities is achieved. With these considerations we obtain the da^'"" subtraction term 
as 



1 ~ S,a 

dCT - 

gg^QQqq 



MloNf ^^T^A^i{Pl,P2,Pi,P^\Pb,P&) XI 



Hl/2){\Mt{lQrig,{M)a,Jg,2Q)\^ + \Ml{lQ,jgA.M)grig,^^^^^^ 

+A0(2Q,?5,3,)|A1°(lQ,i5,%;;t,-,(23)Q)|2jf (pi,pr3,P4) 
+al{3g,ig,4q)\Ml{lQ, (3~4),-; ■,iq,Jg,2Q)\''4^\pi,p2,p^4) 
+A0(lQ,?^,%)|Atg((r4)Q,l,;;3„i^,2Q)|Vf (pr4,P2,P3) 

+a!](4,-Jg,3,)|A^g(lQ,jg,t,;;(3l),,2Q)|2jf)(pi,P2,p54) 

-Gliig,3gAq){\M'i{lQ,mgXj^,2Q)\'' + \Ml,{lQX,m^^ 

-^i?3^lQ,3,,4,-)|A^0((ll)Q,^,,i,,(3lK,2Q)|Vf (;n3,P2,^^4) 
-^£;°(2Q,3,,4,-)|7W^(lQ,i,,j,, (31)^,(2^3)^)127^(^1,^3,^1-4) 



~^0(1q, 3,) (|M°((13)q, i„ 4g- ■ ig, 2q)\' 



MMUmQArr,h'j9,'iQ)\'')4^\m3,P2,P4) 

4^3(2Q,%,4,-)(|>tO(lQ,l,;;3„j„(2~4)Q)p 



+ \MUlQ,'j,,ig-,-,3g,{24)Q)\'')4^\p,,p^3,Ps) 



+^^(lQ,ig,4,-)(^|A^^((14)Q,2Q;;3,,j3,g|^ + |A^^((14)Q,j,,2Q;;3,,g 
-2|7W°((r4)Q,2Q,3,,l,,j-^)|2-(l/2)|Al0((r4)Q,l,;;35,i,,2Q)|2 
-(l/2)|A^0((r4)Q,i3,l,;;3„2Q)p)jf (pr4,P2,^?3) 



+^!J(2Q,i3,3,-)(^|A^^(lQ,(23)Q;;i,-,i3,4,-)|^ + |A^i:(lQ,j,,(23)Q;;i,-,4,-)|^ 
-2|A^0(lQ,(2l)Q,l,-,45,i^)|2-(l/2)|A^0((lQ,4,^;;l,-,j-^,(2'3)Q)p 
-(l/2)|7Wg((lQ,j„4,-;;i,-,(2"3)Q)|2)jf (pi,p^3,P4) 



+ 



iV2 



^£;3°(lQ,3„4,-)|A4g((13)Q,z^,j,„(34)^,2Q)|2jf (pr3,P2,pi'4) 



+^S3"(2Q,3„4,-)|A^^(lQ,i^,i^,(34)^,(23)Q)|2jf (Pi,p53,pi4) 



^;^(1q, i^, 4,-) ( |A^^((14)q, 2q, ; 3„ j„ i,)|^ + |A^^((14)q, j,, 2^; ; 3„ i,) 
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-2|Al0((14)Q,2Q,3„ig,j^)|2)jf)(pr4,P2,P3) 
~AUlQ^g,S,)[\Mlm)Q,2Q;;igJg,4g)\'' + \M^^^ 

-2\Ml{{r3)Q,2Q,tgAg,h)\^)4''\m,P2,PA) 



-2\Ml{lQ,{22.)Q,%A<i,h)?) 4'^ im , P2 , pa) 



~Al{2Q:igA^)(\Ml{lQ,{24)QV,3g39X)^^ 

-2|M°(lQ,(24)Q,3„l„j^)p)jf)(pi,p^4,P3) 



|2 



(5.24) 

Finally, the da^''^ terms are constructed in the usual fashion. They capture the double 
unresolved behaviour of the real matrix element squared when two colour-unconnected 
unresolved partons are present and they remove the double counting of double unresolved 
limits in do"^'". They are given by 

^^It^QQm = -^LoNf (^) §|^d$4(pi,P2,P3,P4;P5,P6)x 

ii,j)eP(5,6) ^ 

~Al{lQ:ig,4g)Al{2Q/jg,3g)\Mlm)Q,{2^^^^ 
-Ul{lQ;ig,3g)A',i2Qjg,Ag)\Ml{ilS)Q,{2A)^^^^^ 



1 



J^3(lQ>^.>4,-)A^(2Q,i,,3,)|A^^((14)Q,(23)Q,J5-,i5)|V^^^(^4,P^3) 



+-^;^(lQ,^g,3,)A;j(2Q,i,,4,-)|A^^((13)Q,(24)Qj,-J,)|V^^^(pr3,pr4) 

(5.25) 



6. Numerical results 

To verify how well the subtraction terms approximate the double real contributions related 
to the partonic process gg — )• QQqq, we have used RAMBO [96] to generate phase space 
points in the vicinity of the singular regions and computed the ratio 

i?=^^p^ (6.1) 

for each of these points. As before, da^^j^Q stands for the double real radiation contri- 
butions while daffj^i^Q is the corresponding subtraction term. In each unresolved limit we 
define a control variable x that allows us to vary the proximity of the phase space points 
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to the singularity. For the difference da^^i^Q — dcj^^^Q to be finite and numericahy 
integrable in four dimensions, the ratio R should approach unity as we get close to any 
singularity. The phase space points were generated with a fixed centre-of-mass energy of 
^/s = 1000 GeV, the heavy fermions were given a mass of 174.3 GeV, and the two hard 
jets were required to have pT > 50 GeV. 

6.1 Double soft limit 

The double soft phase space configurations are characterised by the QQ pair taking nearly 
the full center of mass energy of the event s, leaving the massless final state qq pair with 
almost zero energy, as depicted in fig. 1(a). In fig. 1(b) we show the ratio between the 
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Figure 1: (a) llustration of a double soft qq event, (b) Distribution of R for 10000 double soft 
phase space points. 

double real radiation matrix element and the subtraction term for three different values of 
the control variable x defined in this case by, x = {s — si2)/s. It can be seen that as the 
soft qq pair, takes a smaller share of the total energy, i.e. as x becomes smaller, the peak of 
the distribution around i? = 1 is sharper. This is a sign that the approximation improves 
as the limit is approached. 

6.2 Triple collinear limit 

Other double unresolved configurations in which the amplitude for this process is singular 
are the triple collinear limits ig||3q||4q illustrated in fig. 2(a). In fig. 2(b) we show how, as we 
make the control variable x = — S345/S smaller, that is, as we get closer in phase space to the 
singularity 5g||3g||4g of the real radiation matrix element squared, the histogram becomes 
more pronouncedly peaked around unity. This signals us again that the approximation 
is good. Similar results which are not shown are obtained for the triple collinear limit 
involving the other initial state gluon 6g. 

6.3 Double collinear limits 

The last type of double unresolved limits that the amplitude for the process gg — ?■ QQqq 
has, are the initial-final double collinear limits igW'iq + jg||4q. In fig. 3(a) we represent 
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Figure 2: (a) Ilustration of a triple collinear event, (b) Distribution of R for 10000 triple collinear 
phase space points. 
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Figure 3: (a) Ilustration of a double collinear event, (b) Distribution of R for 10000 double 
collinear phase space points. 



schematically the kinematical configuration where 5p||3g + 6g||4g. As it can be seen in 
fig. 3(b), the convergence of the subtraction term to the real radiation matrix elements is 
indeed achieved as we make the control variable smaller. The same results are obtained 
for the double collinear limit 6^1 |3g + 5g| |4q. 

6.4 Final-final single collinear limit 

In addition to the double unresolved limits discussed above, the double real radiation 
matrix element contains two types of single collinear limits. For the final-final collinear 
limit 3g||4g depicted in fig. 4(a) we obtain the results shown in fig. 4(b), where, once more, a 
good convergence of the subtraction terms to the real radiation contributions is achieved as 
the limit is approached. It should be noted that due to the presence of angular correlation 
terms this good convergence is only achieved after integration over the azimuthal angle of 
the collinear pair. The procedure by which these angular correlations are integrated out 
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Figure 4: (a) Ilustration of a final-final single collinear event, (b) Distribution of R for fOOOO 
final-final single collinear phase space points. 

has been explained in the case where massive partons are present in [69]. 
6.5 Initial-final single collinear limits 

The last type of unresolved limits are the four collinear limits between the massless final 
state (anti) quark and one of the initial state gluons. As it can be seen in fig. 5(b), also in 




(a) (b) 

Figure 5: (a) Ilustration of an initial-final single collinear event, (b) Distribution of R for 10000 
initial-final single collinear phase space points. 

these limits our subtraction terms constitute a valid approximation of the real radiation 
differential cross section. As it is indicated in fig. 5(a) the histogram in fig. 5(b) corresponds 
to the limit 5g||3g. Similar results are also obtained for the other three initial-final single 
collinear limits: 6g||3g, 5g||4g, and 6g||4g. 

7. Conclusions 

In this paper, we present the double real radiation contributions to the tt hadronic pro- 
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duction cross section coming from the partonic process gg — t- ttqq. These contributions 
develop infrared divergencies when one or two partons become unresolved (soft or collinear) 
such that a systematic subtraction procedure is required before these contributions can be 
evaluated numerically. 

We follow the formalism developed in [69], where we extended the NNLO antenna 
subtraction formalism to include the evaluation of hadronic observables involving a massive 
pair of fcrmions. Section 2 contains a summary of the method used in this paper. 

To capture the double unresolved singular behaviour of the gluon-gluon initiated ma- 
trix elements, we derive the appropriate four-parton antenna functions and establish their 
limiting behaviour respectively in sections 3 and 4. Using these, in section 5, we explic- 
itly construct the antenna subtraction terms for the gg — > tiqq subprocess in leading and 
subleading colour contributions. 

In section 6, we check numerically that our subtraction terms approximate the real ma- 
trix elements in all single and double unresolved configurations in a point-by point manner. 
In the regions of phase space associated to the single unresolved collinear kinematical con- 
figurations, the ratio between real matrix elements and subtraction term approaches unity, 
provided the azimuthal terms associated with these collinear limits are correctly treated. 

The double real corrections stemming from the partonic channel gg — )■ ttqq and con- 
tributing to the hadronic production of a top-antitop pair presented here provide a substan- 
tial step towards the calculation of the NNLO corrections to the top-quark pair production 
at the LHC. Future steps include in particular the computation of the remaining double 
real subtraction terms related to partonic channels involving only gluons in initial and fi- 
nal state and the computation of mixed real- virtual contributions for all partonic channels 
involved. 
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A. Three-peirton antennae 

In this appendix we list the three-parton antenna functions used in the subtraction terms of 
section 5 together with their single unresolved limits. The massive and massless antennae 
can be found in [78, 80] and [67, 72] respectively, together with their integrated forms. All 
corresponding single unresolved factors can be found in section 4.1. The massive A-type 
flavour violating initial-final antenna initiated by a gluon, is new. It will be given below in 
unintegrated and integrated form together with its infrared limits. 

In order to make the mass-dependence in the expressions of the antenna functions 
explicit, we use the same convention as in the paper and define our invariants as Sij = 2pi-pj. 
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A.l Massive final- final antennae 

Our subtraction terms use A and E-type massive final-final antennae. The former is given 

by 

^g(lQ,3„2^) = - 1 ^^^ + ^ + ^ + ^ + £13 

(£;2^ + 2m|j \S13S23 Sl3 S23 Sl3 ^23 

^ 2 f 8su 2si2 2si2 2S23 2 2 2si3 

^ V'Sl3S23 Sfg S13 S23 §23 

V -^23 '^13// 
with E^j^ = {pi + P2 +P3)^- It has only a single soft limit: 

A0(lQ,3g,2Q) Si32{mQ,mQ). (A.2) 

This antenna has been derived and integrated in [80] . 
The E-type antenna is 

El{lQ,3g,4q) = -2 f^ia + si4 + ^ + ^ - 2E^rnmQ] + 0{e), (A.3) 

{el. -ml) V ^34 .34 J 

where E'^^ = {pi + P3 + Pa)^ ■ Its only infrared limit is 

^0(1^, 3,,%) J-P^^^a{z). (A.4) 

S34 

This antenna has been derived and integrated in [78] . 

A.2 Massless initial-final antennae 

We use the following A-type antenna 

al{l„ 3„ 2,) = , \ , f-^ + + 2su - ^ - .23) + 0{e), (A5) 

S123 (Sl3 + S23) V «13 Sl3 Sl3 / 

with S123 = S12 — S13 — S23. It is obtained by partial fractionning the full antenna 
A^{lq,3g,2q) SO that it only contains the following single collinear limit 

aO(l„3„25) J-P^^^q(z). (A.6) 

Sl3 

We also need the following G-type antenna 

G'sii,, 3„ 4,) = ^fil + iA+ 0{e) (A7) 

S134 \S34 534/ 

which has the single collinear limit 

3g, 4g) J-Pg^^a{z). (A.8) 

S34 

All massless initial-final antennae have been presented in unintegrated and in integrated 
forms in [72]. 
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A. 3 Massive initial-final antennae 
Initial-Final E and D-type antennae 

The massive initial-final three-parton D-type antenna is given by 

DlilQ,3„l,) = ^--2 fs.is - 9S14 - 6.34 + IfM + fii - fii + ^£HfM 

fQ2^^2y\ S^s .13 .14 ^13 



3.13.34 4.23 2.f4 4.f3 4.f4 6.13.14 2.f3 
H 1 1 1 r Imnm^ 

SlA Sl4 S13S34 S34 S34 .34 .14S34 

, 2 ( 2.^3 4.34.13 4.13 2s|4 6.34 2.I4 4.14 

+"^Q I 2 ' 2 ' 2 ' ' 

\ Sl4 Sl4 Sl4 .14 Sl4 S14S13 .13 

^ 6.34 _ 2sf4 _ 2£|4 _ 4.14.34 _ _ 2m|m^.34 ^ 2mS34 \ ^ 

Sl3 Sl3 Sl3 Sl3 ) S13S14 SlS^M / 

(A.9) 

where = —{p\ + Pz — P'if' and = \fQ^. It has a soft gluon limit as well as a single 
collinear limit 

3g, 4,) ^ 5i34(mQ, 0), (A.IO) 

i?3°(lQ,3„4,) -^P,p^g(^). (A.ll) 
S34 

This antenna is related to the split antennae D\{A,g\ 3^, 1q) and D\(A:g\ 1q, 3g) defined in [78] 
through 

^3(lQ,3„4,) = i?0(4,;3„lQ)-i?°(4,;lQ,3,). (A.12) 



The massive initial-final E-type antenna used in our subtraction terms reads 

ElilQ, 3„ 4,) = i -2 f -.14 + .13 - ^ - — - 2mQm^) + 0(e), (A.13) 

(Q^ + ml) V "34 S34 J 

with = —(pi + P3 — Pa)"^, and = yfQ^- The only infrared limit that this antenna 
has is 

£;30(1q,3„4,) J-P^^^q(^). (A.14) 

S34 

Both of these antennae have been derived and integrated in [78]. 
Massive flavour violating initial-final A-type antenna 

The new three-parton flavour violating A-type antenna initiated by a gluon is denoted by 
3^,25) and reads 

.On 6 9 N 1 / 2.^2 2.12 2.12 , Si3 .23 

^3(lQ>3g,2g) = -Y 1 \ 

[Q^^m^o) V«13S23 Sl3 S23 S23 S13 
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-VHf 



2 / 2.S12 2S23 



'13 



«13 ■513 



+0(6), 



(A.15) 



where = -(pi +p2 -P3)^- 
Its infrared limit is 



^3(l0'3g,2q 



%l|39 1 
S23 



(A.16) 



Its integrated form denoted by ^3(1^,3^,2^) is obtained by integrating the expression 
given in eq.(A.15) over the unresolved intial-final antenna phase space involving one massive 
parton of mass mg [78]. The integrated antenna is given by, 

(Q2 + m2,) 



•^3(lQi3p,2g) 



1 



d$2 



C{e) J 27r 
where C(e), is a normalisation factor given by 



-^3(lQ,3„2,-) 



C(6) = (47r)^ 



87r2 



and the initial-final massive antenna phase space denoted by d$Xj ^j, is given by, 



d$2 



27r 



(A.17) 



(A.18) 



(A.19) 



In eq.(A.19), d$2 is the massive 2 — >■ 2 phase space involving one massive final state 
particle, which we parametrise as [78] 



2 Nl-2e 



2r(l-e)- 

where y runs from to 1. In this context, the invariants are given by, 

Q2 + m2 



2^3, • 

2pi • Pit 



^2 I ^2 



)] 



(A.20) 

(A.21) 
(A.22) 



and the center of mass energy Ecm is 



Er.!!!, — 



'Q^l-x) + ml 



with 



X = 



X 



(A.23) 



The integrated antenna reads 



2pj • q 

1 — 4x + 2x^ + x'X'o + 2x^x0 — 2x^X0 



1 



+ (1 - 2x + 2x^) { ln(l -x) + ^ ln(x^(l - xq)) 



2(1 - xxo) 
+ 0{e) 



with 



Xo 



(A.24) 
(A.25) 
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and the splitting kernel p'qg {x) given by 



p^^^{x) = l-2x + 2x'^. (A.26) 

Having a single collinear limit in its unintegrated form, this integrated antenna develops a 
pole proportionnal to the splitting kernel p~qg {x), as expected. 

A. 4 Massless initial-initial antennae 

Our subtraction terms employ two types of three-parton initial-initial antennae: F and 
G-type antennae. The F-typc antenna is given by 

F3^i„ 2„ 3,) = ^ f - 12512 + 12513 - 12.23 + ^ + 1^ - 

Sl23 V ■513-523 ^13 S23 

_^6S23£12 ^ 6gi3Si2 _^ 4£|3 _ As\^ _^ 25^3 _ As\^ 
Sl3 S23 Si3 S23 SI3S12 S12 

-^-^ + ^-^ + ^\+0{e) (A.27) 
S12 S12 S23S12 / 

with S123 = S13 — S12 — S23- Its infrared limits are 

F°(l^,2^,3^) ^5i23(0,0), (A.28) 

Flilg, 2g, —Pgg^ciz), (A.29) 

F^ig, 2g, 3g) J-Pgg^a{z). (A.30) 
S23 



Finally the initial-initial G-type antenna function is 



1 



GO(i„ 3„ 4,) = — - ^ + 0{e). (A.31) 

with, S134 = si4 — S13 — S34. It only has a single initial-final collinear limit 

G0(i„3„4,) J-Pg^^Q{z). (A.32) 

S34 

Both of these antennae have been derived and presented in their unintegrated and inte- 
grated forms in [72]. 



B. Four-parton antennae 

In this section, we list the two known four-parton antenna functions used to construct 

our subtraction term presented in section 5: the massive final- final B type antenna with 
a massive QQ pair as radiators and the massless initial-initial G-type antenna with two 
initial state gluons as radiators. Their infrared limits were given in section 4. 
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B.l Massive final-final B-type antennae 

This antenna, denoted by ^^(Iq, 4^, 3^, 2q), reads 



-64(1(3, 4g, 3q, 2q) = ^- ^ I \r— [S12S13 + S12S14 + S13S23 + S14S24] 



H 2~ [■Sl2«23 + •S12S24 + •S13S23 + S14S24] 

S34S234 

H 2 2 [2S12S13S14 + S13S14S24 + S13S14S23 - S13S24 - S14S23] 

^34*134 

+ 22 [2S12S23S24 + S13S23S24 + S14S23S24 - S13S24 - ■S14S23] 



H ^- [2si2 + S12S23 + S12S24 + S12S13 + S12S14] 

S34S134S234 

+ -0 [ - ■S13S24 - •S14S23 - •Sl3«24 - S?4S23 

S34S134S234 

+S13S14S23 + S13S14S24 + S13S23S24 + S14S23S24 
-2si2Si3S24 - 2S12S14S23 + 



S134S234 

8gl3gl4 _^ 8523^24 4 ^ 

'*34'*134 '^34'^234 '^134 '^234 



2 2 

2~ [*12 + S23 + S24] 2~ [*12 + *13 + ^u] 

S34Sf34 S34S234 

2 

H [4si2 - Sl3 - - S23 - S24] 

S34'Sl34'S234 

2 ^ [S14S23 + S13S24] 

S34S134S234 



4 



+ 



LS34S134 S34S234 



+ 0(e), (B.l) 



where E'^^ = (pi + P2 + +P4)^- It is normalised to the tree-level two-parton matrix 
element (with couplings and colour factors omitted) 

\M\in* ^ Q0)|' = 4 [(1 - e)El, + 2m%\ . (B.2) 

B.2 Massless initial-initial G-type antenna 

This antenna denoted by ^^(Ig, 3q, 4g, 2g) reads 

G\(^g,\Aq-,%) = [2S13S14S23S24 " s\.^s\^ - s\^s\^ + {s\^ + 8^4] 

^ I '^12*34 ■'I2S13S34 

H [S23 + 2S24S34 + s|4 + s\A + ''^^'^ [Sl4 + 2S24 + S34] 

S12S13S134 S12S13 



_I_£13 1 



S12 S12S34S134 L 

2 3 3 

+^23*24 + ^23 + ^24 



2si4s|3 + 23145^4 + 2S14S23 + 2S14S24 + 523*24 



1 

+ 



S12S34 



2si3Si4 + 4si3S23 + 3si3S24 + 2Si3 
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-SUS23 + 4Si4 
Sl2 



+ 



S12S134 



[2si4S24 - 4S23S24 - S23SU - S24S34] 



S34S134S234 
S12 



[-2si2Si4 + S12 + 4si4S24 + 4Si4 + 4s 



24 



S134S234 



[-6S14 — 6S24 + 3534] + 



2512^14524 

534S134S234 ' 2S13S24 



+ 



SI2S34 



-S14S23] + 



■S24 



[2si2S24 + S12 + 2S24] 



S13S34S234 
— 2si2S24 — S12 + 2523^24 " ^S^^ — 2s\^ 

-s\2 - 2523524 - - S^^ 



S13S34 



2S12S23 



S34 



S13S134 



2si2S23 + 2si2S24 



6S12S24S34 + 4si2S24 



S13S134S234 

+3si2si4 + 3s?2^24 - 3Si2*34 " ■S12 + 3s24si4 " 4s|4S34 + 2^24 - 5^4 



8S12S23 + 4si2S24 - 4S12S34 + 4Si2 + 2523*34 + 6S23 + 25^4 



S13S134 



+25^4 



S13S234 



[4S12S24 - 2S12S34 + S12 - 3S24S34 + 4s|4 + 



3- [2si2 - 2S23 + 2S24 - 2S34] 2^ 

Sl3 S34S134 L 



4si2Sl4 - 8S12S24 + 2Si2 



+4si4S23 + 4si4S24 — 4S23S24 + 4s: 



24 



+ 



'14 



„2 „2 I 
*34*134 '- 



4Sl2S23 + 4S12S24 



-2si2 - 4S23S24 - 2si3 " 2s|4 



+ 



'34 



+4si3S23 - 2S13S24 " 2Si3 + 2si4S23 



4si2Si3 - 4si2Sl4 - S12 
1 



+ 



S34S134 



+5Si2 — 8S14S23 + 6S14S24 + 6Si4 + 6S23 + 4s 



'24 



4S12S23 - 2S12S24 



+ — [2S12 - 2S13 

S34 



1 



-6S14] + -2— [2S12S23 + 2S12S24 - S?2 " 2523*24 " ^23 " •^24] 
*134 

1 7 

+ [4si2 - 4si4 + 4s23 + 2S24 + 3s34] - - + (lo2,3o4) 

S134 2 

+0(e), 

(B.3) 

where = S12 + S34 - S13 - sm - S23 - S24, Si34 = S34 - S13 - sm, and S234 = ^34 - ^23 - S24- 
This antenna is normahsed to the matrix element squared associated to the process gg ^ H 
which is given by 

\Ml{gg^H)\^ = \{l-e)Q\ (B.4) 
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